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TECHNICAL  ACCOMPLISHMENTS 


Introduction 

This  paper  provides  a  summary  of  the  procedures  and  recent  results 
of  a  project  whose  fundamental  goal  is  the  experimental  study  of  unsteady 
turbulent  boundary  layers  in  adverse  pressure  gradients.  Such  boundary 
layers  are  important  in  many  aeronautical  applications,  including  gas 
turbines  and  compressors,  helicopter  rotors,  and  maneuvering  aircraft. 
Presently  a  two-dimensional,  low  Mach  number  flow  over  a  NACA-0012  airfoil 
section  provides  the  test  conditions.  A  rotating  elliptical  cylinder 
located  behind  and  beneath  the  airfoil  trailing  edge  is  the  source  of  an 
unsteady  perturbation  of  variable  frequency.  At  present  operating  speeds 
of  9  to  30  mps,  perturbation  reduced  frequency  ranges  from  0.5  to  6.3. 
Two-dimensionality  of  the  experiment  is  preserved  by  sealing  gaps  between 
the  airfoil  and  sidewalls,  and  by  spreading  the  sidewalls  to  compensate 
for  boundary  layer  growth  on  the  sidewalls  and  maintain  uniform  static 
pressure. 

The  advantage  of  this  technique  is  that  the  airfoil  surface  remains 
stationary,  allowing  hot  wire  velocity  measurements  in  the  boundary  layer 
to  be  readily  obtained.  The  disadvantage  is  that  the  disturbance  amplitude 
is  fixed  for  a  given  frequency,  and  ellipse  location,  and  is  not  uniform 
along  the  airfoil  chord,  being  of  much  greater  amplitude  at  the  trailing 
edge. 

Previous  experimentation  in  this  project  involved  measuring  steady 
state  parameters.  Kanevsky^"  measured  the  pressure  distribution  produced 
by  circular  cylinders  of  varied  diameter  and  position  relative  to  the  airfoil 
and  compared  this  to  the  pressure  disturbance  created  by  the  elliptical 
cylinder  and  to  that  expected  for  a  cambered  airfoil  at  angle  of  attack. 
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Cervisi  measured  pressure  and  boundary  layer  velocity  profiles  for 
different  steady  orientations  of  the  elliptical  cylinder.  This  defined 
the  quasi-steady  behavior  of  the  wing-cylinder  system. 

Induced  Upwash  Measurements 

Research  over  the  past  year  has  concentrated  on  taking  initial  unsteady 
measurements.  This  research  is  reported  in  References  3  and  4.  First 
the  upwash  induced  along  the  airfoil  chord  by  the  rotating  elliptical  cylinder 
was  measured.  An  X-hot  wire  probe  was  used  to  determine  vertical  and 
horizontal  velocity  components.  Figures  1  and  2  show  the  test  configuration 
and  its  relation  to  the  tunnel  cross  section.  The  wind  tunnel  being  used 
is  the  Wright  Brothers  Wind  Tunnel  which  has  an  8x10  foot  elliptical 
cross  section. 

The  instrumentation  used  for  these  initial  tests  was  all  analog  in 
nature,  as  seen  in  Fig,  3.  Each  hot  wire  was  operated  using  constant 
temperature  anemometry  circuits.  The  outputs  were  linearized  and  passed 
through  a  sum/difference  amplifier  to  obtain  vertical  and  horizontal 
velocity  components.  A  trigger  signal  was  produced  by  a  photo-electric 
pulser  on  each  revolution  of  the  elliptical  cylinder  shaft.  This  trigger 
was  used  to  synchronize  real  time  oscilloscope  displays,  FM  data  tape  re¬ 
cording,  and  the  waveform  eductor.  This  device  performed  phase-locked 
averaging  by  dividing  the  periodic  waveform  into  100  segments,  and  time¬ 
averaging  each  segment.  Thus,  for  example,  if  the  cylinder  period  was  AT, 
the  20th  step  in  the  average  was  an  average  of  the  input  voltage  at  time 
equal  20/100  times  AT  after  the  trigger  pulse  over  a  selected  number  of  cycles, 
typically  of  order  100.  The  phase-locked  average  tends  to  eliminate  all 
random  noise  and  turbulence  components  of  frequency  not  equal  to  a  multiple 
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of  the  cylinder  frequency.  This  process  was  primarily  useful  for  determining 


the  amplitude  and  phase  of  the  primary  harmonic  as  well  as  showing  the 
shape  of  the  average  waveform. 

Figure  4  shows  the  mean  induced  angle  plotted  against  distance  from 

the  position  of  airfoil  leading  edge,  for  a  free  stream  velocity  of  30  mps 

coc 

and  reduced  frequency  K  =  ^  varying  between  1  and  5.  Note  the  increased 

OO 

angle  near  the  trailing  edge.  The  increase  with  frequency  is  a  result  of 
increasing  cylinder  circulation.  Figure  5  shows  the  amplitude  of  the 
fundamental  harmonic  against  •  -tance  at  reduced  frequency  equal  to  1  for 
varied  air  speeds.  This  plot  demonstrates  that,  while  minor  variations 
exist  in  induced  angle  of  attack,  reduced  frequency  seems  to  be  the  most 
important  parameter. 

Airfoil  Surface  Pressures 

The  second  series  of  measurements  were  of  the  unsteady  pressures  on 
the  airfoil  surface.  The  upper  and  lower  surfaces  were  each  fitted  with 
17  pressure  taps  along  the  center  line.  The  taps  for  each  surface  were 
connected  to  a  scanivalve  tap  selector  and  capacitive  pressure  transducer. 

By  comparing  the  response  to  a  pure  acoustic  tone  of  this  tap,  tubing,  j  i  ir 

("-i  "  ■  ;  •••<■■  ia 

yarn,  scanivalve,  and  transducer  system  to  that  of  the  transducer  alone, 
the  frequency  response  was  determined  to  be  flat  up  to  800  Hz  with  a  drop 
in  amplitude  to  50%  at  1000  Hz.  Since  the  fundamental  frequencies  considered 
here  range  from  5  to  100  Hz  this  is  felt  to  be  satisfactory.  Figure  6  shows 
the  tap  locations. 

The  instrumentation  block  diagram  for  these  measurements  and  for  the 
velocity  profile  measurements  is  shown  in  Figure  7.  With  the  exception  of 
the  addition  of  a  spectrum  analyzer  and  averager,  the  system  is  the  same 
as  that  discussed  above.  Test  conditions  studied  included  air  speeds  of  9, 
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20,  and  30  mps  with  reduced  frequencies  of  0.5,  1.0,  1.5,  2.0,  4.0,  and 
6.3  .  Initially  the  airfoil  angle  of  attack  was  kept  at  0°.  The  first 
class  of  results  involve  difference  pressures.  Measured  difference 
pressures  were  compared  to  predictions  made  using  Theodorsen's  unsteady 
thin  airfoil  theory’’  based  on  the  previously  measured  induced  upwash. 

Figure  8  illustrates  a  typical  comparison  for  the  mean  pressure  difference. 
The  line  is  the  induced  upwash  based  prediction,  while  the  X's  are  the 
pressure  data.  Figure  9  shows  the  fluctuating  difference  pressure  amplitude. 
Qualitative  agreement  exists,  but  the  measured  amplitude  is  less  than  the 
upwash  prediction.  This  difference  increases  with  reduced  frequency. 

Figure  10  shows  the  phase  shift  of  the  fundamental  frequency.  Agreement 
is  fair  over  the  front  90%  of  the  airfoil,  while  there  is  a  sharp  difference 
at  the  trailing  edge.  Phase  shift  is  defined  here  as  the  difference  in 
time  between  the  minimum  of  the  quantity  and  the  minimum  obstruction 
position  of  the  elliptical  cylinder.  Several  possible  explanations  for 
these  discrepancies  exist.  The  most  important  beinq  the  limitation  of 
the  theory  to  linear,  incompressible  processes  with  the  product  K  M  <<  1. 

In  other  words,  the  effect  of  the  boundary  layer  on  the  unsteady  flow  is 
not  included  in  the  theory.  The  presence  of  a  finite  thickness,  uncusped 

trailing  edge  may  also  contribute  to  the  discrepancy.  In  this  connection 
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Archibald  ,  Satyanarayana  and  Davis  ,  and  Fleeter  ,  have  also  observed 

differences  between  measured  unsteady  difference  pressures  and  predictions. 

The  second  class  of  pressure  results  involve  individual  surface 

pressures.  Mean  surface  pressures,  as  shown  in  Fig.  11,  seemed  quite  well 

behaved,  largely  dependent  only  on  reduced  frequency,  and  then  exhibiting 

only  an  increasing  difference  due  to  the  increased  cylinder  circulation. 

The  fluctuating  surface  pressures  show  more  interesting  behavior.  Figure  12 

shows  the  amplitude  and  Fig.  13  the  phase  for  the  fundamental  at  reduced 
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frequency  0.5  .  These  results  are  qualitatively  similar  to  the  behavior 
expected  based  on  the  quasi-steady  results.  One  point  of  interest  is  the 
large  jump  in  phase  at  the  leading  edge  stagnation  point.  A  major  shift 
in  behavior  occurs  for  reduced  frequencies  of  2  or  larger.  Figures  14  and 
15  show  amplitude  and  phase  of  the  fundamental  at  reduced  frequency  4.0  . 

The  amplitude  is  somewhat  reduced,  particularly  at  the  leading  edge,  but 
the  primary  difference  is  in  the  phase  over  the  front  20%  of  the  airfoil. 
Instead  of  the  gradual  increase  in  upper  surface  phase  seen  at  lower 
frequencies,  a  sharp  drop  occurs.  Large  phase  shifts  near  the  front  on 

9 

unsteady  airfoils  were  also  seen  by  Franke  and  Henderson  ,  No  convincing 
explanation  of  this  behavior  has  yet  been  found. 

Two  summary  plots  are  shown  in  Fig.  16  which  shows  the  behavior  of  the 
fluctuating  pressure  coefficient  on  the  upper  and  lower  surface  and  their 
difference  at  X/C  =  0.7.  The  important  point  to  note  is  the  two  distinct 
zones.  At  low  reduced  frequencies  one  type  of  behavior  occurs,  namely  the 
pressure  coefficient  is  constant  or  decreasing,  while  at  high  reduced 
frequencies  the  pressure  coefficient  is  increasing.  This  implies  that  there 
is  a  change  in  the  dominant  phenomena  from  the  convective  processes  at  low 
reduced  frequencies  to  smaller  scale  local  unsteady  processes  at  higher 
reduced  frequencies.  In  contrast  to  this  result.  Fig.  17  shows  that  the 
RMS  pressure  coefficient  change  across  the  airfloi  tends  to  approach  zero 
at  the  trailing  edge.  The  nature  of  the  approach  seems  only  weakly 
dependent  upon  the  reduced  frequency.  Table  I  shows  the  relation  of  these 
results  to  other  results  in  the  literature. 
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Airfoil  NACA  64A010  NACA  65  Circular  Arc  NACA  0012  NACA  64A010  NACA  0012 
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SUMMARY  OF  UNSTEADY  AIRFOIL  DATA 


Boundary  Layer  Velocity  Profiles 


The  third  series  of  unsteady  measurements  were  airfoil  boundary  layer 
velocity  profiles.  The  physical  set  up  is  seen  in  Fig.  18.  The  hot  wire 
probe  may  be  translated  perpendicular  to  the  surface  in  0.02  mm  or  greater 
increments,  with  a  total  travel  of  9  cm.  Translation  is  performed  using 
a  DC  motor,  with  relative  probe  position  being  determined  with  a  linear 
potentiometer.  Since  the  boundary  layer  thickness  (99%  of  external  mean 
velocity)  is  typically  6  mm  at  70%  of  chord  and  1.25  cm  at  94%,  the 
6  cm  travel  allows  flexibility  in  positioning  the  probe  at  various 
chordwise  locations  and  airfoil  angles  of  attack,  while  keeping  probe  motion 
normal  to  the  surface.  Some  problems  have  occur ed  in  operating  within  the 
nearest  0.2  mm  from  the  surface.  Vibrations  caused  by  both  wind  tunnel 
fluctuations  and  by  the  rotating  cylinder  can  cause  the  hot  wire  to  brush 
against  the  surface  and  break  at  the  joint  to  its  support.  The  inability 
to  reliably  approach  the  surface  has  also  affected  the  determination  of 
the  position  of  zero  height.  It  had  been  hoped  to  find  zero  height  by  the 
completion  of  a  low  current  circuit  when  the  probe  touched  the  surface, 
however  not  being  able  to  safely  get  near  the  surface  prevented  this.  Zero 
height  was  operationally  determined  by  extending  the  linear  velocity 
profile  in  the  laminar  sublayer  to  zero  velocity.  Heights  determined  in 
this  manner  were  consistent  to  within  0.06  cm  for  varied  profiles  at  the 
same  location. 

Figure  19  shows  a  mean  velocity  profile  for  K  -  2.0,  X/C=.94  .  Mean 
profiles  tended  to  be  relatively  insensitive  to  the  unsteady  fluctuations. 
That  is,  the  mean  profile  for  K=1  is  very  similar  to  that  for  K=6.3  for 
the  same  chordwise  location  (and  therefore  about  the  same  mean  pressure 
gradient) . 
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Figure  20  shows  a  logarithmic  mean  velocity  profile.  Note  that  the 


axes  are  referenced  to  the  external  velocity  rather  than  to  the  friction 
velocity.  This  is  because  methods  to  estimate  the  friction  velocity 
from  the  velocity  profile  have  not  given  convincing  and  unambigous  results. 
An  unsteady  skin  friction  gauge  is  being  developed  to  independently  find 
the  friction  velocity,  but  it  was  not  available  for  these  tests.  The 
roultiplicitive  factor  in  the  velocity  axis  is  a  streamline  curvature 
correction.  This  correction  is  the  cause  of  the  downturn  in  the  profile 
for  logarithmic  heights  above  4,  which  actually  are  above  the  boundary 
layer. 

Figure  21  shows  a  root-mean-square  velocity  profile.  The  RMS  velocity 
includes  the  fundamental  frequency,  its  harmonics,  turbulence,  and  wind 
tunnel  fan  noise.  Figure  22  shows  the  fundamental  harmonic  amplitude. 

Note  in  particular  the  factor  of  2  increases  in  amplitude  perturbation 
velocity.  External  amplitudes  ranged  from  5%  of  the  mean  velocity  at 
95%  of  chord,  reduced  frequency  1.0  or  less  to  0.5%  for  the  68%  of  chord, 
K=6.4  case.  The  velc^ity  phase  shift  relative  to  the  external  velocity 

is  shown  in  Fig.  23.  Lower  frequency  cases  in  general  had  smaller  shifts, 

of  order  5  or  10  degrees,  while  the  high  frequency  cases  were 
characterized  by  ^Positive  shifts  of  15°  or  higher  in  the  logarithmic 
region  of  the  mean  profile,  switching  to  negative  phase  near  the  surface, 
in  the  region  of  highest  velocity  amplitude. 

Future  Work  v. 

v 

While  the  above  results,  produced  largely  using  analog  instrumentation 
are  encouraging  in  that  they  define  the  basic  characteristics  of  the  wing- 
rotating  cylinder  system,  future  work  including  the  collection  of  data 
for  non-zero  airfoil  angles  of  attack  and  for  additional  chordwise  boundary 


layer  measurement  locations  will  require  the  implimentation  of  a  digital 
data  collection  and  real-time  analysis  system.  This  system  will  involve 
triggering  the  data  acquisition  using  the  cylinder  shaft  pulser,  digitizing 
and  storage  of  unprocessed  data,  phase-locked  averaging  of  data  over  many 
cycles  (similar  to  that  previously  done  using  the  eductor)  Fourier 
transformation  to  find  the  amplitude  and  phase  for  arbitrary  harmonics, 
ensemble  averaging,  and  real-time  graphics  display  of  selected  results 
to  verify  equipment  operation  and  to  quickly  identify  new  phenomena 
deserving  further  study. 

Also  to  be  added,  as  mentioned  above,  is  a  wall  shear  gauge.  This 
will  be  important  both  to  normalize  properly  the  velocity  profile  and 
because  wall  skin  friction  is  in  itself  an  important  aerodynamic  parameter. 

Future  near  term  research  will  involve  the  study  of  unsteady  pressure 
and  boundary  layer  profiles  at  non-zero  angle  of  attack,  up  to  static 
stall  and  velocity  profiles  across  the  wake  to  attempt  to  shed  light  onto 
the  unusual  behavior  of  the  phase  of  the  trailing  edge  difference 
pressure.  Further  away  are  studies  of  unsteady  separation,  its 
correlation  to  other  unsteady  boundary  layer  processes  including  the 


evolution  of  discrete  structure  within  the  boundary  layer. 
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SPOKEN  PAPER  AT  64TH  SEMI-ANNUAL  MEETING 
SUPERSONIC  TUNNEL  ASSOCIATION 
9-10  October  1980 

NEW  DISCOVERIES 

The  results  generated  in  this  contract  period  include: 

1.  While  surface  pressures  at  the  trailing  edge  are  difficult  to  measure, 
data  up  to  X/C  =  0.98  extrapolates  such  that  the  Kutta  condition  seems 
to  be  satisfied  up  to  reduced  frequencies  of  6.4  .  (Note  we  define 
the  Kutta  condition  as  "The  velocity  remains  finite  at  the  trailing 
edge".  Thus  a  finite  jump  in  pressure  at  the  trailing  edge  satisfies 
the  Kutta  condition  under  this  definition) .  The  phase  lag  at  the 
trailing  edge  is  geater  than  predicted  by  simple  theory. 

2.  In  the  Reynolds  Number  range  of  .3  to  1.0  x  10^  and  in  the  MIT 
Wright  Brothers  Wind  Tunnel  the  unsteady  pressure  distribution  seems 
to  be  well  correlated  as  a  function  of  reduced  frequency. 

3.  Aerodynamic  Interference  is  suitable  for  producing  an  unsteady  forcing 
function  that  will  produce  a  well  behaved  unsteady  boundary  layer, 

as  supported  by  preliminary  data. 
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DIMENSIONS  IN  CM 


WIND  TUNNEL  TEST  SECTION 
FOR  INDUCED  UPWASH  TEST 


1 .  ROTATING  ELLIPTICAL  CYLINDER 

2.  X  HOT  WIRE  PROBE 

3.  PITOT  STATIC  PROBE 
H.  2  HP  DC  MOTOR 

5.  FLYWHEEL 

6.  PHOTOELECTRIC  PHASE  SENSOR 

7.  SIDEWALLS  AND  SUPPORTS 


FIG.  2 
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PRESSURE  AND  VELOCITY  TEST  IK  S  TRUMSN  TAT  I  OR 
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CP  DELTA  CP  LO  CP  UP 


FLUCTUATING  PRESSURE  AMPLITUDE 
U  INF  -  20  MPS  X/C  =  0.70 
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CP  RMS  A  CP  RMS 


RMS  DIFFERENCE  PRESSURE 
NEAR  AIRFOIL  TRAILING  EDGE 


UINF  CMPS)  :  +03,  *20,°30 
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CONFIGURATION  FOR  PRESSURE 
AND  VELOCITY  PROFILE  STUDIES 


V13H1/A 

FIG.  19 


■PUT  * 


33 


UMEAN/UE;  UINF=  ZO  K=  2.0  X/C=.94 


jgjD  □□□q 


LOGCY+UE/NID;  UINF»  20  K=*  6.4  X/O.94 


CVi  (S3  CD  CD  <-  OJ  QQ 

'ri  "H 


VJ.3H1/A 

FIG.  2  / 


35 


URMS/UE  UINF=  20  K=  6.4  X/C=94 


in 

rvi 


<s 

rvi 


in 


S> 

•H 


in 

Q 


s 

<si 


V13H1/A 

FIG.  22 


36 


UFLUCT /UEFLUCT ;UINF“  20  K<=  2.0  X/O.S4 


